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We consider an exactly solvable model for production of fermions in the Friedman flat universe 
with a scale factor linearly growing with time. Exact solution expressed through the special functions 
admit an analytical calculation of the number density of created particles. We also discuss in general 
the role of the phenomenon of the cosmological particle production in the history of universe. 
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I. INTRODUCTION 

Although a cosmological particle production [T-|4| is 
well known phenomenon p-Qj exactly solvable models 
are of interest because they would allow compar- 
ing the different approximate calculation methods for the 
number of produced particles. One of the problems in 
this field is to define the in- and out- vacuum states if 
the metric has not an asymptotic form corresponding to 
that of the Minkowski space-time [U [ll[ . 

As the approximate methods, the adiabatic vacuum 
[j|, WKB series [13] and instantaneous Hamiltonian di- 
agonalization Q methods are widely used. Evidently, 
the later method does not define an exact vacuum state. 
As is shown in Ref . fi~3| , a certain superposition of the 
states, which diagonalize the instantaneous Hamiltonian, 
is needed to obtain the exact vacuum state. 

The method has been suggested in Ref. |l4| , which al- 
lows finding the exact vacuum state by minimization of 
some functional. Here, on basis of this method, we have 
found the exact vacuum states for the case of a flat uni- 
verse with a linearly growing scale factor. Such a model 
has appeared in Ref. (TBj . where the universe driven by 
the vacuum has been considered. Although, that model 
in its present version contradicts the nucleogenesis the- 
ory, which insists on an existence of the radiation epoch, 
it is of theoretical interest as a toy model. It should be 
noted, that the scalar particle creation in such a Universe 
was considered earlier [l6j (See e.g. Ref. [13]) on basis of 
the WKB approximation, and father exact solution have 
been obtained [lfj . For fermions energy-momentum ten- 
sor have been calculated in Ref. 181. 
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Let us also emphasize, that the flat universe with a 
linearly growing scale factor differs substantially from the 
Milne-like empty closed universe in a sense that the later 
can be transformed to the region of the Minkowski space- 
time (l9l ] and does not exhibit a particle production. 



II. FERMION PRODUCTION 

Let us consider the fermion in the expanding universe. 
After decomposition of the bispinor ip(r,rj) in the com- 
plete set of modes ij)(r) = ^ fe ipkirf) e lfer , the Lagrangian 
of the fermion field in the expanding universe 20] takes 
the form 



k 

-a 3 ip£(ak)if; k - a A rmp£f3ip k , (1) 

where r\ is the conformal time, a(rj) is the universe scale 
factor. 

The equation of motion is 
Set' 

Wk - { ak )i>h + ' l ~^-i>k - ma/3ipk = 0, (2) 
2a 



The fermion field is quantized as 



i>h = btt. s v-k,s + a k , s Uk, s , 



(3) 



where the bispinor is: 



iXfc + maxk / <Ps 

. Ma® ips 
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Uk,s(lj) 



where ip s are ip+ = \ q 1 ano - V- 
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The bispinor Vk, s is expressed as Vk, s — i7°7 2 (ufe,s) T , 
where the symbol T denotes a transposed vector, u — 
?i + 7° and the representation of the Dirac matrices is the 
same as in Refs. [U ^ 
M 



The functions Xk(v) satisfy 



where u±(k), v±(k), U±(k), V±(k) denote u±.k, u±,fe, 
W±,fej V±,fc with k directed along z-axes. From Eqns. (O 
it follows that the functions Xk, %k corresponding to the 
different vacuums are connected as 



Xk + (fc 2 + m 2 a 2 - ima') Xfe = 0, 
fc 2 XfcXfc + (amxt - ix'k*) {amXk + ix'k) = 1> 

that corresponds to a time-dependent oscillator with a 
complex frequency. 

Linearly growing scale factor 



X k = cosr fc Xk(v) 
K = cosr k x.Uv) - 



-iSi 



smr k {ma{rj)xt{rj) ~ iXkiv)) / k > 
sinr fe (ma(i])xk(v) + ix'kiv)) A- 



According to [14|, the quantity 



<rk(v) = 7;(x' k (v)x*k(v) + Xl' (v)Xk(v)) 



(6) 



a(t) = a H t 



(4) 



corresponds to the exponential dependence a(rj) — 
aoexp("Hr/) in the conformal time dr/ — dt/a, where the 
conformal Hubble constant is T~L = a^H^,. To describe the 
particles production, one should define the in- and out- 
vacuum states. There is a lot of solutions connected by 
the Bogoliubov transformation (23| . 

For the fermions with the momentum k directed along 
z-axes, the Bogoliubov transformations are: 

U+{k) = cosrfc — sinrfce* 5 '" v—(— k), 

V+(fc) = cosr-fe V+(k) + sinr k e~ lSk w_(-fc), 

U-{k) = cosrfc-u_(fc) — sinrfce l5fc k), 
V-(fc) = cosrfc«_(fc) + smrke~ lSk u + (—k), 



(5) 



is monotonic in the past but has oscillating behavior in 
the future for the in-vacuum. For the out-vacuum state, 
it is oscillating in the past but is monotonic in the future. 

With the help of the functional minimization method 
[3], we obtain that 



/csch(27r/c/'H) , H , 

Xk = V^M —i ex P ( _jm a o e /rl - ikr)) 



Uk 



L^fi?ima Q e^/n) (7) 



for the in-vacuum state and function a shown in Fig. [T] 
(a) has the form 



Ok = 



""-< *»aech (iff] ( Ji * ( enVma o) Jlk 



iWt 



H 



H 



I H 2 



e^mao 
H 



J _ £fe _ 1 

n 2 



For the out-vacuum state (see Fig. [T] (b)) 

X k = — exp (nk/2H - ima Q e Hr) /H - ikr]) 
( ik-H 2ik 2ie nv ma 



(9) 



Here U (a, b, c) is the confluent hypergeometric Kummer's 
function and L^' (z) is the generalized Laguerre function 



Square of sin r*fc can be expressed as 



sin 2 Tk 



k { Xk X' k - Xkx'k) 



im 2 a 2 XkX* k - ma X* k x'k + ma XkX* k ' + ik2 XkX* k + i-XkXt 

I 



1 + exp(27rfc/"H)' 



(10) 



The density of the produced fermion is 



III. DISCUSSION AND CONCLUSION 



N 



Va 3 



sin 2 r k 



Vd 3 k 



3 W 
40-7T 5 a 3 
3 



40tt ; 



H*(t). (11) 



It is interesting to compare the number of the pro- 
duced fermions with the result of Refs. [f| [25j], where the 
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FIG. 1: (a) Dependence of the function <7fc on the conformal time for the in- vacuum state and (b) for the out- vacuum state; 
m = 1, k = 2, H = 2. 



dependence 



a(t) = a t q = a rf, 
p = q/(l-q), 0<q<l 



(12) 



(13) 



was considered. The density of the produced fermions as 
well as scalar particles have been estimated as @, [25| 



n ~ m 3 (mty 3q . 



(14) 



Under q = 1 one has n ~ f -3 . Up to a numerical 
multiplier this coincides with our exact result given by 
dU, (|TU|) . (ITT1) . The numerical multiplier equals approxi- 
mately to ^5 « 2.5 x 1CT 4 . 

Now we want to discuss the moment of time when the 
particles of some mass m are created. According to Ref. 
[14| , the moment of time, from which the function a given 
by (J6j> acquires oscillating behavior corresponds to the 
particle creation. 

From the analysis of the expression (JSJ) it follows that 
the particles of mass m are created by time 771 deter- 
mined by the equation moo exp (Hr/i) /% ~ 1 that is 
m/H{t\) ps 1. By that moment of time according to 
(fTTj) total density of a particles approximately equals 
ni ~ H(ti) 3 (or equally n\ ~ m 3 ) and further it is 



simply reduces with time as the n(t) 



From 



dimension arguments it can be argued that for arbitrary 
a(t) density number of the created particles can be esti- 
mated as 



n(t) 



,3 a 3 (h) 



a 3 (t) 



(15) 



Let us check this formula for the dependence given by 
(IT21 . For the time of particle creation we have ti « — , 
and finally come to (fl"4")l . 

It is also interesting to discuss in general the possible 
role of the cosmological particle production in the cre- 
ation of matter in the Universe. According to the modern 
view decelerating fast-roll stage preceded the inflationary 
stage |26| . At this stage Hubble constant decreased down 
to the value H(t\) right up to the time t\ when the in- 
flation begins. Let us estimate amount of matter created 



before the inflation stage. At the inflation stage typi- 
cal value of the Hubble constant is Hi — m, where m 
is the inflanton mass. It is an order of m ~ 10 _6 M p , 
where M p ~ 10 18 GeV is the Plank mass. Density of 
the inflantons produced under the background of the co- 
herent inflanton field is give by (|15p . During inflation 
and further expansion their density is reduced in a 3 /a 3 , 
times where ao is the scale factor of the present universe. 
Present time matter density created due to this effect can 
be estimated as 



mn\- 



(M p 10" 
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(16) 



where for simplicity we suggest that the created inflan- 
tons decay to the massive particles, so that their gen- 
eral amount of mass is conserved. Comparing this quan- 
tity with the observed density of matter in the universe 
p c = Mp Hq = Mp 1CT 122 we find that the expansion rate 
has to be 2a. = 10 33 in order to the densities be equal. 

That is the matter arising due to cosmological creation 
of inflantons will be negligible only if the expansion rate 
of the universe during inflation and further expansion 
exceeds 10 33 ( 76 e-foldings). 

To summarize, we have presented the new exactly solv- 
able model for fermion production in the flat universe 
with a scale factor linearly growing with time. Spec- 
trum of the produced fermions over physical momentums 

p = k/a 



exp(2vrp/iJ(t)) + 1 



is thermal at large p, and the Hubble constant divided by 
2-7T plays a role of the temperature. The number density 
of the created particles does not depend on the particle 
mass, and particles of large mass can be easy created. 
However, the function <Jk begins to oscillate later in the 
time, when the mass decreases. 

The oscillations of this function mean that the real 
particles appear [lij]. Thus, in the limit of m — > 0, the 
particles would be created infinitely late. 
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